Scaling of the elastic contribution to the surface free energy of a nematic on a 

sawtoothed substrate 
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We characterize tlie elastic contribution to the surface free energy of a nematic in presence of a 
sawtooth substrate. Our findings are based on numerical minimization of the Landau-de Gennes 
model and analytical calculations on the Frank-Oseen theory. The nucleation of disclination lines 
(characterized by non-half-integer winding numbers) in the wedges and apexes of the substrate 
induces a leading order proportional to q\nq to the elastic contribution to the surface free energy 
density, q being the wavenumber associated with the substrate periodicity. 
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I. INTRODUCTION 

Topological defects are ubiquitous in many branches of 
physics, spanning from condensed matter physics [l|, Q 
to cosmology 0] . They may emerge in systems with bro- 
ken continuous symmetry, and their presence can intro- 
duce essential singularities in the free energy which lead 
to infinite-order phase transitions such as the Kosterlitz- 
Thouless transition [1,[H. Liquid crystalline phases, such 
as nematics, are prototypical for showing topological de- 
fects, such as point hedgehogs and disclination lines [EH- 
Topological arguments show that within the bulk and in 
3D the only stable disclination lines have ±1/2 winding 
number The presence of substrates may induce the 
nucleation of disclinations with other winding numbers. 
Theoretically, non-half-integer disclination lines are pre- 
dicted to nucleate in the cusps (wedges and/or apexes) of 
substrates [sl-ITlj , as well as in sharp boundaries between 
domains characterised by different anchoring conditions 
in flat surfaces [H, Furthermore, there is some ex- 
perimental evidence of the formation of these unusal non- 
half-integer disclination lines on surfaces [l^ . 

It is well known that the nematic director field in pres- 
ence of structured substrates may be distorted, leading to 
an elastic contribution to the free energy. Since the sem- 
inal work by Berreman [l^ this problem has been 
extensively studied and generalized in the literature 0- 
[Tll. [l7l - [2l] |. However, most of these studies focussed on 
smooth substrates or on the effect of isolated cusps. In 
this paper we will consider the effect that the disclina- 
tion lines nucleating on the cups of the substrates have 
on the elastic contribution to the surface free energy of 
the nematic on periodic and cusped structures. In par- 
ticular, we will consider a sawtoothed substrate which fa- 
vors homeotropic anchoring. Other geometries may also 
be considered with our formalism, such as step-like sub- 
strates, which have been studied numerically and exper- 



imentally in the context of zenithal bistable switching in 
nematic devices [13, 112, HI] ■ We will show that the discli- 
nation lines induce a contribution to the elastic contribu- 
tion to the free energy density (per unit projected area) 
which depends only on the geometric characteristics of 
the substrate, and that scales with the wavenumber q 
associated to the periodicity of the substrate as —qlnq. 

The paper is organized as follows. The different models 
for the nematic in the presence of a sawtooth substrate 
arc presented in Section II. We focus in particular on 
substrates which favor homeotropic anchoring and their 
consequences. In Section III we report the numerical 
and analytical results obtained from these approaches. 
In Section IV we present our conclusions. 



II. THE MODEL 

We consider a nematic phase in contact with a saw- 
tooth substrate characterized by an angle a and a side 
length L (see Fig. [T]). At the substrate the nematic 
molecules preferentially align homeotropically, i.e. par- 
allel to the local normal to the substrate. The system 
is translationally invariant along the out-of-plane axis z 
and periodic along the x axis. We will only consider az- 
imuthal distortions for the nematic director n, which can 
be parametrized in terms of the angle 9 between the di- 
rector and the y direction as n = (— sin 6, cos 9, 0). Out- 
of-plane or twist deformations may also be important un- 
der other conditions, as a twist instability may occur [23 |. 
but we checked numerically they are not relevant for our 
choice of parameters as free-energy minimization always 
leads to azimuthal distorted textures. Finally we impose 
that, far from the substrate, the bulk nematic phase ori- 
ents homogeneously either along x (the iV" texture) or 
the y direction (the texture), which are the only rele- 
vant situations allowed by symmetry considerations. We 
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study the orientational ordering of this system within two 
different models: the Landau-de Gennes (LdG) model 
and the Frank-Oseen (FO) model. 

A. The Landau-de Gennes model 

In the LdG model, both isotropic and nematic phases 
can be locally represented by a traceless, symmetric 
order-parameter tensor with components Qij , which can 
be represented as Qij = ^S[ninj — ^Sij] + ^B[lilj—mimj], 
where are the Cartesian components of the director 
field n, S is the nematic order parameter which mea- 
sures the orientational ordering along the nematic di- 
rector, and B is the biaxiality parameter, which mea- 
sures the ordering of the molecules on the orientations 
perpendicular to n, characterized by the eigenvectors 
1 and m. The LdG free energy can be written as 

J^hdG = /v(0bulk + 0el)dy -I- /_4(Asurf ds wllCrC 0bulk IS 

the bulk free energy density, (/)ci is the elastic free energy 
density, and (psuri is the surface free energy are defined 
as [11]: 

,/,buik = aTrQ2-feTrQ3 + c[TrQ2]2 (1) 

0oi = ^dkQijdkQij + ^djQijdkQik (2) 
2 

<f>surf = ^ 3^ Tr[Q • Qsurf] (3) 

where a depends linearly on the temperature, b and c arc 
positive constants, and Li and L2 are positive parameters 
related to the elastic constants. If we will rescale all 
the variables as follows [2^: Q = 6cQ/6, the positions 
f = r/^, where the correlation length ^ is defined as 
e - 8c(3ii + 2L2)/b^ and AdG = 24^03 J■LdG/C'&^ we 
get that AdG = /v(0buik + 4>ci)dV + J^4>sur{ds, with 
rescaled free energy densities: 

0buik = ^TTrQ2_^TrQ3 + l[TrQ2]2 (4) 
<^oi = [dkQijdkQij + KdjQijBkQtk] (5) 

6 + ZK 

2 ~ ~ 

0surf = --wTrlQ- Qsurf] (6) 

Here t = 24ac/6^ is a dimensionless temperature, 
K = L2/L1 is an elastic dimensionless parameter (k > 
—3/2) and w — 16wc/b^£, is the dimensionless anchoring 
strength. Hereafter we will consider these rescaled ex- 
pressions, so we will drop the tilde notation. For r = 1, 
the bulk free-energy density has two minima correspond- 
ing to (/>buik = for rescaled scalar order parameters 
S*! = (isotropic phase) and S'n = 1 (nematic phase), so 
both phases are at coexistence. It is important to note 
that the order parameter 5* in the coexisting nematic 
phase is rescaled, so its value in real units is &/6c, which 
must be smaller than 1 (typically « 0.4). If the elastic 
parameter k is positive (negative), the nematic prefers 



Nematic phase 



y 




FIG. 1: Schematic picture of the geometry of the system, 
characterized by the side length L and the angle a. 

to align parallel (perpendicular) to a possible nematic- 
isotropic interface. Finally, Qsmt defines the favored ten- 
sor at the substrate. We will favor a homeotropic align- 
ment of the nematic by setting Qsmf = (3i>' (E) — l)/2, 
with v the normal vector to the substrate, establishing a 
direct connection to previous papers |26l - [29j . Hereafter 
we will restrict ourselves to the nematic phase at the 
nematic-isotropic (NI) transition (r = 1), with k ~ 2. 

B. The Frank-Oseen model 

The FO model can be considered as an approximation 
to the LdG model, in which we assume that the varia- 
tions of the nematic order parameter arc restricted to the 
neighbourhood of the substrate (with a width of order of 
^) and inside defect cores, and takes the bulk value S 
elsewhere. Assuming there are no disclination lines in 
the bulk, substitution of this ansatz into the LdG free- 
energy functional leads to the following FO approach for 
large L (in units of ^): 

^ « J-po = 0buik(5)V + / dsJ:{0{s)) + ^ /dy|V0p 

Ja ^ Jv 

(7) 

where V is the total volume of the nematic, A is the 
substrate area in contact with the nematic, the elastic 
constant K = (9/2)52(2 + k)/(3 + 2k) and T,{e) is an 
effective anchoring potential to be determined in the fol- 
lowing that is due to the nematic order parameter dis- 
torsions close to the substrates. 

In order to make a quantitative comparison with the 
results within the LdG model, instead of considering an 
ad hoc expression for the anchoring potential, we will de- 
rive it from the LdG model. As it was mentioned above, 
the nematic order parameter distortions are assumed to 
be confined to a layer of width 77 ~ ^ close to the sub- 
strates. On the other hand, the variations of the nematic 
director field close to the walls are rather small in di- 
rections parallel to the surfaces for large L (except close 
to the wedges and apexes). So, we may estimate S(6') 
as the LdG excess free energy per unit area of a slab of 
width 7], for which, at the boundary y = 0, we impose 



a surface field (/)suif, and for y > r] we consider a bulk 
nematic phase with an uniform nematic director char- 
acterized by a tilt angle 6 with respect to the y axis. 
Fig. [5] shows obtained by numerical minimization 

with a conjugate-gradient method for rj = 1.5^ (simi- 
lar results are obtained for other values of 77). For each 
value of w, the minimum value of corresponds to the 
homeotropic alignment = 0, where it takes the (true) 
nematic- wall surface tension value (7nw(w')- For not too 
large values of 9, T,{6) takes a Rapini-Papoular form [s^l 
T,{6) ^ (Tnw + ^oiw) sin^ 9. There are deviations for 9 
around tt/2, but in any case the involved energies are 
much larger. In the range of values of w we will consider 
in this paper (0 < w < 2), both ctnw and Sq are of or- 
der of w. Indeed (tnw can be obtained analytically from 
minimization of the LdG functional in the presence of a 
flat wall and for homeotropic anchoring [2^ . In this 
situation the equilibrium nematic director field does not 
show any deformation and there is no biaxiality, so the 
resulting free energy at NI coexistence can be expressed 
in terms of the nematic order parameter profile S = S{y) 
as .Flcig = Aj^iS^ - + S^ + {S'f/2)dy - wS{Q). 
where S'{y) = dS/Ay. By introducing a magnetization- 
like field m{y) = 2S{y) — 1, the LdG functional reduces to 
the Landau-Ginzburg free-energy functional of an Ising 
model for zero applied magnetic field in presence of a flat 
wall. This problem has been extensively studied in the 
literature jsil [s^ , leading to a nematic order parameter 
profile S{y) = (1 + g exp(— v^y))"^, where g is obtained 
from the boundary condition S"(0) = —w. The resulting 
expression for ctnw is: 



(Tnw 



\/3(g + 3)g2 w 



6(1 + 3)3 1 + g 



(8) 



where g = -{1 + l/y/2w) + \J {I + 1/V2w)2 - 1. 

In order to obtain the azimuthal angle 9 field, we have 
to minimize the energy functional given by Eq. ([7]). In 
the bulk , we thus have to solve the Laplace equation 
V'^^ = with appropriate anchoring conditions at the 
boundaries. Let v = (— sin6'o, cos6'o, 0) now be the lo- 
cal normal to the substrate. Using the approximation 
E(0) ^ (Tnw + Eg(w)sin^(0 — ^o) found previously, the 
nematic director field then satisfies approximately that 
isTiv- V6'-|-So sin(2(6'-6'o)) = 0. Note that we may assume 
strong homeotropic anchoring conditions when L is large 
compared to the extrapolation length K/2Yiq ^ K/w, 
or equivalently wL » 1. We can justify this assump- 
tion by considering the following rescaling: r* = r/L, 
9*{v*) ~ 9{r) and J"pQ = Fyo/L. In this rescaled de- 
scription, we must minimize the FO free energy ([T]) in the 
rescaled domain subject to an effective anchoring poten- 
tial Y,*{9*) = LY,{9). This means that the dependence 
on the size of the system can be absorbed into the co- 
efficients of the anchoring potential: (Tnw (which does 
not affect the anchoring conditions) and Sq (which are 
rescaled by a factor of L). Strong anchoring condition 
is then satisfied when the rescaled extrapolation length 
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FIG. 2: Plot of the effective anchoring potential S(6) as a 
function of sin^ 6, for a layer width 77 = 1.5^ and w — 0.5 
(circles), w = 1 (squares), w = L5 (diamonds) and w = 2 
(triangles). The dashed horizontal lines correspond to the 
values of ctnw obtained from Eq. and the continuous 

lines correspond to the fit to a Rapini-Papoular expression 
for small values of the tilt angle. The values of the slopes Eq 
are 0.5655 {w = 0.5), 0.8729 {w = 1), 1.0678 {w = 1-5) and 
1.177ii {w = 2). 



K/2T,qL ^ 1, in agreement with our previous estimate. 



III. RESULTS 

A. Numerical results within the LdG model 

We now turn to the evaluation of the distortion contri- 
bution in the LdG model We consider values of a between 
and 7r/2, 8 < L < 96 and < w < 2. In order to get 
the orientational ordering, we numerically minimize the 
LdG free-energy functional by using a conjugate-gradient 
method. The numerical discretization of the continuum 
problem is performed with a finite element method com- 
bined with adaptive meshing in order to resolve the dif- 
ferent length scales that may emerge in the problem [s^ . 
The numerical procedure is similar to that used to study 
the wetting transition by the nematic phase in this geom- 
etry psl ] . The numerical minimization of the LdG func- 
tional J^^dG shows that the N-^ texture has lower free 
energy if a < 7r/4 owing to lesser distortion. Conversely, 
the iV'l texture has lower free energy for a > 7r/4. These 
results are in agreement with previous studies in the liter- 
ature 

As a consequence, we observe bistability in 
a range of values of a around 7r/4. Fig. [3] displays some 
typical textures. We see that on the substrates the ne- 
matic orientation is preferentially homeotropic. Inspec- 
tion of these textures show that in general there are no 
disclinations in bulk, except for large w and a around 
7r/4, but there is an important distortion on the nematic 
director field close to the wedges and apexes. On the 
other hand, when disclinations are observed in bulk (see, 
for example. Fig. [3] (b) or (c)), they are at a distance of 




FIG. 3: (Color online) Contour maps of the nematic order 
parameter S (blue (dark grey) for small S, white for the bulk 
value S = 1, and orange (light grey) for higher S) and the 
nematic director field n (grey segments) for the equilibrium 
textures obtained from minimization of the LdG model for 
L = 16 and w = 1.0: (a) a = tt/6 (N^), (b) a = tv/4 
(N^), (c) a = 7r/4 (iV") and (d) a = 7r/3 (Af"). The singular 
nematic director field from the FO approach (Eqs. (|ll|l and 
(|12|l ') are also shown for comparison (black segments). Note 
that they are indistinguishible almost everywhere. 

order of ^ from a wedge or apex, which is almost inde- 
pendent of L. 

The analysis of the calculated equilibrium free ener- 
gies confirms that the leading order contribution to the 
equilibrium free energy per unit area is ctnWj which is 
a further confirmation that there is strong anchoring on 
the substrates. As a consequence, we may obtain the 
elastic part of the free energy per cell (i.e. a slice paral- 
lel to the yz plane with width along the x axis equal to 
the substrate period) and unit length in the z direction, 

/elastic; '^S 

2L(JxdG|eq/-^ — cnw)- Numerical evaluation 
of /clastic shows a clear dependence on L in a broad range 
of values of w (see Fig|4]for a = tt/6 and a = tt/3). It is 
interesting to compare this result to the prediction within 
Berreman's approach for the elastic contribution to a 
smooth substrate- nematic surface free energy [lB|. Al- 
though originally this approach was introduced for shal- 
low sinusoidal substrates with strong tangential anchor- 
ing, it can be extended straigthforwardly to the present 
case with homeotropic anchoring, either weak or strong. 
By using the electrostatic analogy, we may expand the 
azimuthal angle field 6 as: 

oo 

0{x, y) = ^ M„ sm{nqx) cxp{-nqy) (9) 

n=0 

where q is the wavenumber associated to the substrate 
periodicity, and the coefficients M„ are chosen to satisfy 




FIG. 4: (Color online) Plot of /dastic as a function of the an- 
choring strength w for a — n/6 (black symbols) and a — n/S 
(orange (light grey) symbols). We considered different cell 
sizes; L = 8 (circles), L = 16 (squares), L — 32 (diamonds), 
L — 48 (triangle up), L = 64 (triangle left) and L — 96 
(triangle down). 



the imposed boundary conditions. Substitution of this 
expression in the FO free energy leads to the elastic con- 
tribution to the surface free energy which depends only 
on Aq, with A the roughness amplitude. For sinusoidal 
substrates under moderate or strong anchoring condi- 
tions, Berreman's approximation /elastic ~ TTK{Aq)^ is 
valid if qA is smaller than 1 [g^] . In any case, if A scales 
like L, the surface free energy should be independent of L. 
In order to explain and characterize the anomalous scal- 
ing for the elastic free energy obtained within the LdG 
framework, we will now resort to the FO approximation. 



B. Results with the FO model 

Disclination lines are known to nucleate on wedges 
and apexes 0-[ll|. Their presence introduce a singu- 
lar contribution to the orientation field 9 which can- 
not be expanded as Eq. ([9]). We may write 0{x,y) as 
9s{x,y) -\- 9^s{x,y), where 9s is the singular contribution 
due to the disclination lines, and 9ns is the non-singular 
part. Note that Berreman assumed that 9s = 0. The sin- 
gular 9 field associated to one isolated disclination line 
located at the origin and charaterized by a winding num- 
ber / takes the form — / arctan(a;/y) = /5(ln(i(^)), where 
C, ^ X + iy and 3(C) is the imaginary part of Q. In the 
bulk, the values of the winding number / are restricted 
to half-integer values to avoid discontinuities in the ne- 
matic director field. For a periodic linear array of discli- 
nation lines, we may obtain the singular nematic orien- 
tation field from a conformal transformation: 
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where the wavenumber which characterizes the substrate 
periodicity is g = -k / L cosa. The transformation maps 
a neighbourhood of any point C„ — n{2Lcosa), with 
n integer, to a neighbourhood of the origin in the 
complex plane. If we denote C = Cn + £ (kl small), 
then the transformation leads to a complex number 
C w (-l)"ge/2. Thus the orientation field /3(ln(iC')) = 
— /arctan(tan(ga;/2) coth(g?//2)). which is a solution to 
the Laplace equation, periodic in x with period 2Lcosa, 
and singular at each reduces to the field corre- 
sponding to an isolated defect of winding number / in 
the neighbourhood of any C„. Focussing on the range 
X € [—i cosa, L cosa], the azimuthal angle generated 
by this function is = at a: = and x = ±Lcosa. 
On the other hand, for \qy/2\ ^ 1, 6 « — sgn(j/)/qx/2 
for X 7^ ±icosa, where sgn(j/) — y/\y\ is the sign of 
y. Note that there is a jump from 9 ~ — sgn(?/)/7r/2 to 
9 = sgn(y)/7r/2 when crossing x = Lcosa. In the bulk, 
this solution is acceptable if / is integer (note that the 
nematic state is invariant under inversion of the nematic 
director). If / is half-integer, we may add to this solution 
the orientation field / arctan(tan(ga;/2)) due to an ar- 
ray of disclination lines with winding number — / located 
a-t Cn — ioo, which in the x interval (—L cosa, L cos a) 
reduces to Iqx/2. The compound orientation field 
/[— arctan(tan((7x/2) coth((j'j//2)) -I- arctan(tan((7x/2))] 
has no discontinuity at x = ±Lcosa for y > (in fact, 
it goes to zero as qy/2 3> 1), but it jumps from Itt to 
—It: when crossing x — Lcosa for y < 0. This is again 
physically acceptable if / is half-integer. Finally, we must 
mention that there are different physically equivalent rep- 
resentations of the same singular field associated to a pe- 
riodic array of disclination lines located at x = (,„ as, 
for example, the field /[arctan(cot((j'a;/2) tanh((jy/2)) — 
arctan(cot((3'a;/2))]. This orientation field shows a phys- 
ically acceptable discontinuity at 2; = for y < if / is 
half-integer. 




FIG. 5: Plot of the integration contour C (dashed line) to 
get the elastic contribution to the free energy in the region 
R inside the contour. The contour is deformed to avoid the 
cores of the dischnation lines with wandering numbers Ii and 
I2 at the wedge and apex (filled circles). We highlight the 
contour Ci with an associated outwards normal v. 



wavenumber q and for x G (—L cosa, L cos a) as: 
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(12) 
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When disclinations lines are located on surfaces, their 
winding numbers / are not constrained to half-integer 
values [sMllj. as we may tune the winding numbers in 
order to match the boundary conditions close to any 
wedge or the apex in the strong anchoring regime. Alter- 
natively, a Schwartz- Christoffel transformation may be 
used Q. For the N-^ texture, we find that the wind- 
ing numbers Ii and I2 for the disclination lines at the 
wedge bottom and at the apex top, respectively, are 
= -a/(7r/2 - a) and h = a/(7r/2 + a). For the 
A^ll texture, the topological charges are Ii = 1 and 
I2 = — (7r/2 — a)/(7r/2 -I- a). The resulting singular ori- 
entation fields, 9^ and 9\ , can be written in terms of the 



arctan 
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This expression does not show any discontinuity above 
the substrate. We checked that these solution approxi- 
mate quite well the boundary anchoring conditions (but 
not exactly), so that (?ns may be neglected in most of the 
cases. Furthermore, the agreement between our ansatz 
for the director orientation field and the numerical LdG 
textures is excellent except close to the defect cores, as 
can be seen in Fig. [3l 

The elastic energy can be obtained by standard con- 
tour integration techniques. By symmetry, /clastic can 
be obtained as twice the contribution of the half-cell re- 



gion R (sec Fig. [5]) . By using the identity V • {0'V9) = 
I V6'p + 9W^d, and using the divergence theorem and the 
fact that \/^6 = in R, we obtain that /elastic can be 
expressed in terms of an integral over the contour C of 
R: 



elastic 



K (t> Ov- veds 



(13) 



Along the segments x = 0, x = Lcosa and y — > oo, 
we impose that 9 = Osmg ~ otoo, where the far-field 
azimuthal angle is aoo = for the N-^ texture and 
ttoo = 7''/2 for the A^" texture. We link Ci to the x ~ Q 
and X — Loos a segments by circle arcs of small radii e. 
These circles take into account the presence of a defect 
core, and their radii are proportional to the nematic co- 
herence length, which provides a microscopic cutoff for 
the continuum Frank-Oseen model. In these boundaries, 
approaches, up to a constant, the field of one isolated 
disclination line, which satisfies the following property 
u _L V0 as e 0. So we can neglect the contribution 
of the circle arcs to the contour integral. Finally, 9 = a 
along the segment C'l for strong anchoring conditions. 



Taking into account that 
rewrite Eq. ([T^ as; 



TC 



clastic 



K 



K{a 



u ■ V6'ds = 0, we might 



a^) 6 V- VOds (14) 



If 9 is split into a the singular and non-singular contri- 
bution, we see that the non-singular term, which can be 
expressed by Eq. ©, leads to a contribution to /elastic 
independent of L. On the other hand, the leading order 
contribution to /elastic comes from the singular orienta- 
tion field, namely the contribution close to the disclina- 
tion lines and apex. If s is the distance of a point of the 
contour Ci close to the wedge or apex, i'- Vising ~ —Ii/s 
or — /2/s, respectively. The leading order contribution to 
/elastic comcs from integration on Ci as: 



/elastic ~ IC{a) Ini/e 
where IC is defined as: 



lC{a)=K{a^a^){h~h) 



a < 



(15) 



(16) 



Expressions Eqs. p5)) and ([T6|) show that the equilibrium 
texture is N-^ for a < 7r/4, and A^" for larger values of 
a, in agreement with our LdG calculations and previous 
results reported in the literature [Ij. 



C. Analysis of the results within the LdG model 

Finally, in order to check the accuracy of our approx- 
imation, wc analyze the results obtained with the LdG 
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FIG. 6: (Color online) Plot of B{w,a) as a function of the 
anchoring strength w for a = n/6 and a — tt/3. The meaning 
of the symbols is the same as in Fig. 3] 



model by fitting the values of /elastic obtained from min- 
imization to an expression: 



/elastic =/C(a)ln(L/C) + B(a,«;). 



(17) 



where /C(a) is given by the expression Eq. (jl6p and 
B{a,w) is expected not to depend explicitely on L if wL 
is large enough (i.e. under strong anchoring conditions). 
The resulting curves for different values of L of B(a, w) 
as a function of w for a fixed value of a collapse into 
a master curve (see FiglGj. Deviations only appear for 
small values of u; < L~^, at which the strong anchoring 
condition breaks down. 

This result shows clearly that there is a contribution 
to the elastic part of the surface free energy which scales 
logarithmically with the periodicity of the substrate. For 
large L (small q) and w ^ L~^, the surface free energy 
density (i.e. per unit projected area on the xz plane) 
of a nematic in the presence of a sawtooth substrate has 
the asymptotic behaviour ctnw / cos a— {IC{a) /2T:)q lnq + 
0{q). The non-analytical contribution — qlng introduces 
a slow decay of the surface free energy of the nematic 
with increasing L. This may help to explain the large 
deviations with respect to the Wenzel law observed for 
the wetting transition by nematic of a sawtoothed sub- 
strate in contact with the isotropic phase j28j . 



IV. CONCLUSIONS 

In this paper we have analyzed the size dependence of 
the elastic contribution to the interfacial free energy den- 
sity of a nematic in presence of a sawtoothed substrate. 
The nucleation of non-half-integer disclination lines of 
the nematic director field on the apexes and wedges of 
the substrate, predicted in the FO model, induce a non- 
analytical contribution which scales with the substrate 
periodicity wavenumber q as —qhiq for small q. This has 
been confirmed by the numerical results in the full LdG 
model. The periodicity scaling of the elastic contribution 
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to the surface free energy is different from that obtained 
for smooth surfaces, which scales hnearly with q. Our 
arguments are not specific to this kind of substrate, and 
can be extended straightforwardly to any surface shape 
which shows ridges, cusps and similar singularities. On 
the other hand, the size-scaling of the elastic free energy 
has consequences for surface transitions such as wetting. 
Finally, the recent impressive advances in microfluidic 
technology and surface patterning open the posibility of 
an experimental verification of our predictions. 
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